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o ; Abstract 

P^ ' The relaxation of inagnetically confined plasmas in a toroidal geometry is analyzed. From the 

d • equations for the Hermitian inoments, we show how the system relaxes towards the mechanical 

equilibrium. In the space of the parallel generalized frictions, after fast transients, the evolution of 
collisional magnetically confined plasmas is such that the projections of the evolution equations for 



'^ ■ the parallel generalized frictions and the shortest path on the Hermitian moments coincide. For 

O . spatially-extended systems, a similar result is valid for the evolution of the thermodynamic mode 

(i.e., the mode with wave- number k = 0). The expression for the affine connection of the space 

(N 

^ ' covered by the generalized frictions, close to mechanical equilibria, is also obtained. The knowledge 

0\' 

lO . of the components of the affine connection is a fundamental prerequisite for the construction of the 

(N 

^^ \ (nonlinear) closure theory on transport processes. 
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I. INTRODUCTION 

It is known that a macroscopic description of thermodynamic systems requires the formu- 
lation of a theory for the closure relations. A thermodynamical field theory (TFT) has been 
established in order to determine the (non linear) deviations from of the Onsager coefficients 
for thermodynamic systems far from equilibrium [l| . The nonlinear transport equations have 
)een derived by imposing that the thermodynamic theorems for systems out of equilibrium 
2| and the De Donder-Prigogine principle, also refereed to as the Thermodynamic Covari- 



ance Principle (TCP) [3| (see the footnote [l9[), are respected 



Magnetically confined tokamak plasmas are a typical example of thermodynamic systems 
out of Onsager's region. In this case, even in absence of turbulence, the local distribution 
functions of species (electrons and ions), as well as the distribution function for the fiuctu- 
ations of the thermodynamic quantities, deviate from the (local) Maxwellian. According to 
the Onsager theory of fiuctuations, in this condition, the thermodynamic fiuxes are not (in 
general) linearly connected with the conjugate forces (ref. to the Onsager theory [4] and, 
for example, [5]. See also the end of section UTTj) . In tokamak plasmas, the thermodynamic 
forces and the conjugate fiows are the generalized frictions and the Hermitian moments, 



Is 



respectively [6]. The neoclassical theory is a linear transport theory (see, for example, J6|) 
meaning by this, a theory where the moment equations are coupled to the closure rela- 
tions (i.e. fiux-force relations), which have been linearized with respect to the generalized 
frictions (see, for example, Ref. [7|). This approximation is in contrast with the fact that 
the distribution function of the thermodynamic fiuctuations is not a Maxwellian and could 



B^Iq 



provoke some disagreements with the experimental profiles |8[, [9|. It is, however, important 
to mention that it is well accepted that the main reason of this discrepancy is attributed 
to turbulent phenomena existing in tokamak plasmas. Fluctuations in plasmas can become 
unstable and therefore amplified, with their nonlinear interaction, successively leading the 
plasma to a state, which is far away from equilibrium. In this condition, the transport prop- 
erties are supposed to change significantly and to exhibit qualitative features and properties 
that could not be explained by collisional transport processes, e.g. size-scaling with machine 
dimensions and non-local behaviors that clearly point at turbulence spreading etc. (see, for 



example, Ref. |10J ) 



The nonlinear transport equations have been largely used for studying transport processes 



in non equilibrium systems such as magnetically confined plasmas, materials submitted to 
temperature and electric potential gradients or chemical reactions. In magnetically confined 
plasmas, the nonlinear transport equations provide a link between the generalized frictions 
(the thermodynamic forces) and the Hermitian moments (the conjugate flows). This allows 
to determine the particle fluxes (electrons and ions) and energy losses as well as the (non- 



linear) particle distribution 
the papers reported in Ref. 
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iinctions. This task has been accomplished in Ref. [8| and in 
and in the footnote J20l |. 
Finally, the nonlinear transport equations have been derived by adding a further assumption: 
There exists a thermodynamic action, scalar under TCT, which is stationary for general 
variations in the transport coefficients and the affine connection of the thermodynamical 
forces space |lj. However, the determination of this action requires the knowledge of the 
affine connection. The expression of the affine connection can be derived by analyzing several 
examples of relaxation. 



In the manuscript reported in Ref. [ll|, we analyze the relaxation of chemical reactions to 
stable steady-states in the chemical affinities space. In these cases, the small parameter e 
measures the distance of the system from the steady state. 

In this paper, we shall be concerned with the relaxation of magnetically confined plasmas in 
a toroidal geometry. The characteristic feature of the evolution equations is the presence of 
a smallparameter e, the drift parameter, defined as the Larmor radius over a macroscopic 
length p]. In this situation, the long-time behaviour of the solution, describing the evolution 
of the system near the steady state, may be obtained by using the multiple time-scale 



perturbation expansion (see, for example, the book cited in Ref. 12] )• Starting from 
the balance equations for mass, energy and higher order hermitian moments, applied to 
magnetically confined plasmas, we show, in Section[Tll the validity of the following theorems: 

(a) After fast transients, in the Onsager region of the generalized frictions space, a homo- 

geneous system relaxes towards the mechanical equilibrium along a straight line. 

(b) After fast transients, out of the Onsager region, a homogeneous system relaxes towards 

a stable mechanical equilibrium such that J^ U^{Xev.-tr., q) — Jpi i^^i^s.path, q) = 0{e^) 
where W^i^X^v.-tr., q) and U^{Xspath, q) are the evolution equation for the generalized 
frictions and the shortest path equation in the generalized frictions space, respectively. 
X^ and J^ denote the vector of the generalized frictions and the vector of the Hermitian 



moments, respectively. The trajectory traced out by the system and the shortest past 
are parametrized by q. This parameter is defined in the Subsection III Al 

We assume that the mechanical equilibrium is a steady-state of the plasma. Similar theorems 
for the relaxation of the thermodynamic m,ode (i.e., the mode with wave-number k = 0) to 
the mechanical equilibrium, can also be derived for spatially-extended plasmas. 
Even though, the example examined in this paper refers to the relaxation of magnetically 
confined plasmas in a toroidal geometry, the results obtained are valid generally because 
the dynamics include all relevant moment equations and the parameter e is not related to 
the distance of the system from the stationary states. Such example enable us to determine 
the expression of the affine connection for the thermodynamical forces space, near the non- 
equilibrium steady-states. This task is accomplished in Section IIIII 

II. RELAXATION OF MAGNETICALLY CONFINED PLASMAS 

In this section we analyze the case of magnetically confined plasmas. For simplicity, we 
consider fully ionized plasmas defined as a collection of electrons and positively charged ions 
[TJ . The balance moment equations for mass, heat and the non-priviledged fifth hermitian 
moments read, respectively (6| 

T^dtqf'^ = n^T^ermnq'^^Hn + T^Q^'^ + g^^^^ + ^"^'^ + 0{^) (l) 

where cff denote the Hermitian moments of the distribution functions, Qr indicate the 
dimensionless generalized collisional friction terms and gr , (jr are the dimensionless 
source terms. The suffix a distinguishes the two electric species i.e., a = e for electrons and 
a = i for ions. Index n takes the values n = (1, 3, 5). Moreover, Q^ and r^ are the Larmor 
frequency and the relaxation time of species a, respectively, ermn and br denote, respectively, 
the Levi-Civita tensor and the unit magnetic field components. In this paper, the summation 
convention over repeated indices is understood. On the contrary, in Eqs ([T]), as well as in 
what follows, there is no summation over the index a. Eqs ([T]) are the moment equations 
describing collisional plasmas in the presence of an inhomogeneous and curve magnetic field. 
These plasmas are referred to as (magnetically confined) tokamak plasmas. In this paper we 



deal with collisional transport processes (i.e., not turbulent plasmas) where collisions are the 
only source of irreversibility or dissipation. In this case, we can distinguish two collisional 
transport regimes: the fully collisional transport regime (or the Pfirsch-Schliiter transport 
regime) and the low-collisional transport regimes (or the plateau and the banana regimes). 
Physically, these two transport regimes are distinguished by the order of magnitude of 
\mfpi the mean free path. When the mean free path is much smaller than the macroscopic 
length scale Lm we are in the so-called short mean free path regime or, the fully collisional 
transport regime (fctr). The macroscopic length is associated with the gradient lengths of 
the macroscopic quantities, such as the density gradient, the temperature gradient, etc.. Lm 
is defined as the shortest of these lengths. On the other hand, if the fusion temperatures 
are realized, the mean free path A^/p can easily be made to exceed Lm- In this case we are 
in the so-called long mean free path regime or the low-collisional transport regime (Ictr). In 
this regime, the orbits of the particles are bounded in the space (the banana orbits) (for easy 
references, see for example Refs J6|). A useful transport theory must cover both the fully 
collisional and the low-collisional transport regimes. By a magnetized plasma it is meant a 
plasma where the scale length characterizing the plasma (in general the hydrodynamic scale) 
is much larger than the gyro-radii of its constituent charged particles. Thus, by defining 
with pl the Larmor radius (the radius of the circular motion of a charged particles in the 
presence of a uniform magnetic field), the plasma is magnetized if the parameter 

..f. (2) 

is much less than one. Wherever is valid the approximation e ^ 1, we are within the so- 
called drift parameter approximation. It can be shown that plasmas in the both transport 
collisional regimes, described above, can be validly treated in the drift approximation [see, for 



f. B]. 



example, ref. [6|]. The presence of the small parameter e allows to neglect in the equations 
all terms of order e^ or higher. In Ref. |6| an intrinsic kinetic form of the dimensionless 
(density of) entropy production a" of species a is derived under the sole assumption that 
the state of the plasma is not too far from the reference local equilibrium state. It is shown 
that, in the local dynamical triad [see Fig. ([T])], the entropy production is closely associated 




FIG. 1: The local dynamical triad. The radial direction e^ is the unit vector radially directed 
pointing outwards. The parallel direction eii coincides with the unit magnetic vector i.e., eii = 
B/5. The right-handed coordinate system is obtained defining the perp-tangential direction e/\ as 

Ga = ep Aeij. 



with the collision term and it can be brought into the form 
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For easy reference, we list here explicitly the relations between the dimensionless and the 
corresponding dimensional Hermitian moments: 
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^ = (P>II>A) 



where ma, Uq, and T^ are the mass, the number density and the temperature of species a, 
respectively. Moreover, jr, F" and Q" indicate the electric current, the particle fluxes and the 
heat fluxes, respectively. L" is the dimensional fifth-order Hermitian moment corresponding 
to Qr ■ For completeness, we also report the relation between the pressure tensor vr"^ and 
the second-order tensor Hermitian moment qr 
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The dimensionless source terms are defined as 
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with r = (p, II, a) and g^i " are the traceless tensorial Hermitian moments. Indicating with 
e the absolute value of the charge of the electron and with Z the charge number of the ions, 
we have Cq, = — e for electrons and Cq, = +Ze for ions. E and u are the electric field and 
the centre-of-mass velocity of the plasma, respectively. Moreover 6„ = B/i? and ermn is the 
Levi-Civita symbol. Eqs ([3]) allows us to identify the thermodynamic forces X"^ and the 
thermodynamic flows J^: 

X"^ = g;^""^ + -g^^"^^ ; j; = gj^^^^ /i = (1, 2, 3), n = (1, 3, 5) (7) 

Hence, the thermodynamic forces are related to the spatial inhomogeneity and they are ex- 
pressed as gradients of the thermodynamic quantities whereas the thermodynamic flows are 
the parallel component of the quantities defined in Eq. (jlj). Different, but also very impor- 
tant, quantities are the transport fluxes defined as the radial fluxes averaged over a magnetic 
surface. These fluxes are denoted by < qp >, where < ■ ■ ■ > indicates the magnetic- sur- 
face averaging operation. In this notation, < g^ -^ > (i.e. n = 1) and < Qp > (i.e. n = 3) 
denote the (dimensionless), magnetic averaged, radial particles flux and radial heat flux, 
respectively. Let us stress that the transport fluxes are different from the conjugate ther- 
modynamic flows. This main conceptual step deserves some additional explanations. The 
thermodynamic forces, given defined in Eq. (JTj), are not directly accessible to measurement 
in a tokamak. Only surface-averaged quantities, depending solely on the radial coordinate p 
are experimentally accessible. Typically, experiments provide us with temperature, density 



or pressure profiles and with transport fluxes. The characteristic feature of the neoclassical 
expression of the surface-averaged entropy production for toroidally confined fully collisional 
plasmas, is 

< ^e >fctr = < qf^ >fctr ^J'^''+ < qf^ > fctr of^ 

< O, >fctr = < qf^ >fctr gf^ (8) 

where gp is the dimensionless radial pressure gradient (denoted by VpP) 

sr=-.(^)"'— v.i' (9) 

Thus, the average entropy production, depending on the parallel components of the entering 
quantities, is expressed entirely in terms of the transport fluxes and radial forces. This 
inversion of roles is a direct consequence of the zero- divergence geometrical constraint. The 
latter transforms the average entropy production into a radial contribution. However, the 
thermodynamic form shown by Eqs ([H]) is also due to the fact that, in the fully collisional 
regime, there is no contribution of the unprivileged moments, q?} j^ to the average entropy 
production. This is not the case for the low-collisional transport regime. As a result, in this 
latter regime, the average entropy production does not derive from a thermodynamic form. 



which would have been expected [13 1 



< Oe >lctr t <qp >lctr fi'p + < qp' ' >lctr gp" + gyictr 9\\ 

< (J^ >lctr + < qf^ >lctr gf^ (10) 

where ql, and ^m are the average dimensionless parallel electric current density and the 
average dimensionless induced external parallel electric field (denoted with -Ei^), respectively 

-(1) / ^ (1) ^ (1) - (^eTi\^/'^ j(l) e(l) 

9,r=(|^)'\^<BB,f> ; A = ^7b^ (U) 

In literature, this phenomenon is referred to as the divorce between entropy production and 
transport theory for the banana transport mechanism. Further physical and mathematical 



details concerning this subject can be found in the paper cited in the footnote 2l|]. The 



transport fluxes are related to the parallel components of the dimensionless generalized 
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frictions and the dimensionless source terms by [see, for example, Ref. [6|)] 

< qp >fctr= —KaTa < ^f 1 — ^]Q?\ > in the fully colhsional regime 

< Qp >ictr= Ka < j-g'u > in the low-collisional regime 

In general, the quantity Ka is not a surface quantity but it becomes one for some configura- 
tions of the magnetic field as, for example, in the standard model. As previously mentioned, 
Eq. (fT2|) is a direct consequence of the zero- divergence constraint, introduced by the confin- 
ing geometry, which couples the flux components in a new way and introduces the effect of 
the parallel thermodynamic forces on the transport fluxes. The evolution equations of the 
moments, Eqs ([T]), have a hierarchical structure: the equations for the third moments will 
involve the fourth moments, and so. The dynamics of a thermodynamic system is finally 



based on the set of balance equations Eqs ([T]) coupled to the closure equations [22|], relating 
the dimensionless generalized collisional friction terms to the thermodynamic fluxes. In this 
example, the closure equations can be cast into the form (7| 

r.Q"^"^ = -r"^^'"^ (m, n = 1, 3, 5) (13) 

where Qmn (for simplicity, we omit the suffix a) is a positive-definite matrix, which may 
depend on the thermodynamic forces X"'^. In case of magnetically confined plasmas, this 
matrix is symmetric [6] . Inserting Eqs ( !T3l) into Eqs ([T]) we obtain the balance equations for 
the thermodynamic forces 

TadtJ^ = -Sp^r^J: - ^"1 + 0(e2) where (14) 

{0 if /i 7^ K ^ ^ 

(15) 
1 if /i = K 

The task of the transport theory is to describe the evolution of the gradients (of density, 
temperature, magnetic field etc.) towards global mechanical equilibrium. As long as the 
step length in the random walk is smaller than the gradient length scale this process is 
diffusive, so the associated time derivative is expected to be of the order 

(16) 

where D ~ Pi/Ta denotes the diffusion coefficients. Thus, for a magnetized plasma, very 
close to the steady-state, the partial time derivative in Eq. ( 1T4|) is assumed to be small in 
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dt 




J^AT 




Tr, 



the sense 



Q 



-1-^^ (17) 



The weaker version 

r4~e (18) 

describes the relaxation of the plasma, evolving towards the steady-state on a shorter time- 
scale. To sum up our discussion, during the relaxation, close to the steady-state, we require 
that 

i) In the transport regimes (i.e., neglecting fast transients), the time derivatives r^dtJ"^ are 
assumed to be of order e^. During the relaxation, near the steady-states, TadfJ^ are 
of order e. The weaker versions of the rate of the change in time represent the fast 



E 



evolution of the system [6|. 
ii) The Universal Criterion of Evolution theorem (UCE), established by Glansdorff- 



Q. 



Prigogine [2|], and the covariance properties o 
to the thermodynamic forces transformations 



the evolution equations, with respect 
l|, are satisfied. 



Finally, we recall that we are allowed to neglect in the equations all terms of order e^ or 
higher. 

A. Homogeneous Plasmas 

In this case, Eqs (IT4|) may be rewritten, in a more convenient way, as 

(5^ + g^'^X'^'-g^^^^y^g.^X"'' = - 5^.(rV," - X"'') + 0{e') (19) 

where 5^ denoting the Kronecker delta and we have introduced the transport equations re- 
lating the thermodynamic flows with the thermodynamic forces that produce them, through 
the matrix of the transport coefficients, g^^ (we omit the suffix a) : 

J: = g,,X'^^ (20) 

In principle, the matrix of the transport coefficients may depend on X"^. Quite naturally, 
in solving this equation close to the steady-state, we should make an optimum use of the 
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presence of the small parameters. In a magnetized plasma, in accord with experiment, the 
gyrofrequency Q exceeds the collision frequency z/, 

n « 1 (21) 

Taking into account that the mean-free path A = vti^ and pl = vt/^, where Vr indicates 

the thermal velocity, we obtain 

A = ^ (22) 

A 

So, Eq. (12T]) again expresses the smallness of the Larmor radius. A condition, which almost 
always pertains in practice, is 

^ ~ ' P3) 

Physically, Eq. f l23|) implies that, when the plasma is very close to the steady-state [see con- 
dition i)], the time derivative in Eq. ( IT6l) is considerably smaller than the transit frequency. 

^ - ^'— (24) 

Eq. fl23j) is synonymous with u/Q ~ 0{e) and admits as a special case u/Q ^ e. Finally, the 
fundamental expansion parameter is that of the smallness of the Larmor radius pi compared 
with the macroscopic scale length Lm- Hence, the problem boils down to seek a solution for 
X°'^, which may be written as a perturbation expansion in powers of e, 

^a^. = Xo"'' + eX^" + e^X^ + ■ ■ ■ (25) 

where X^^ denotes the reference state. The definition of this state is given below. From 
here on, we are tempted to start with the standard perturbation theory. But this approach 
induces singularities. Indeed, some of the contributions to the correction terms may be 
proportional to the positive powers of time, t". It follows that, in spite of the smallness 
of the parameter, e™ for long enough times, these terms become very large and invalidate 
the classification Eq. fl25p . This fact is so much more troublesome considering that we 
are interested in a solution close to the steady-state i.e., in the long-time behaviour of the 
solution. This difficulty may be overcome by using the multiple time-scale perturbation 
expansion (see, for example, the book cited in Ref. 13]). Thus, we artificially extend the 
unique time variable t into several dimensionless variables g-i, qq, qi, . . . , defined by the 
differential equations 

1 . , 

Ta Q-1 = -, Ta ^0 = 1, Ta Ql = 6, ■ ■ ■ (26) 

11 



Once the multiple time-scale solution has been obtained in this form to the desired degree 
of accuracy, we go back to the original time variable by the replacements 

1 ^ i t 

Q-i^ , Qo^—, Qi^e—,--- (27) 

^ ^CK ^a ^a 

The time derivative becomes 

did d d , . 

" dt e dg-i dgo dgi 

The "singularity" introduced by the term 1/e represents the evolution on the shortest time- 
scale whereas the dependance on ^o represents the slower evolution on a longer time-scale, 
etc. Finally, the solution X"'^(t) is split into many other new variables X"^(^_i, qq, Qi,- ■ ■) 
(with i = 0, 1, 2, ■ ■ ■ ) depending on many (dimensionless) times ^_i, ^O) Qir ' ' ■ As men- 
tioned, the idea is based on the observation that, with the ordering p^ ^ A^/p and pi <^ Lm, 
there exist widely separated time-scales in the system. This allows to define the reference 
state Xq^ as the state depending only on the shortest time-scale typically, the inverse Larmor 
frequency in our problem. Hence, the reference state satisfies the equations dX^^ /dgi = 
(for i = 0, 1, ■ ■ ■ ). The other variables X"^ (with i = 0, 1, ■ ■ ■ ) may, however, depend on a 
longer time-scales. The unknowns X"^ (with i = —1, 0, 1, ■ ■ ■ ) are completely determined by 
imposing the Fredholm solvability (or, the solvability conditions for differential equations) 
on each iteration. This formal procedure allows to eliminate the secular terms in each order 
in e and expansion Eq. fl25l) is uniformly valid for all times because it takes into account 
not only the size of the terms X"^, but also their rate of changes in time. According to 
Eq. flTTl) . on the time-scale ^2? the evolution is so slow that we can consider the system at 
the steady-state, X^l'g-i-^^^, defined by the equations dX^^^gf-g^f-^/dgi = (for i = — 1,0,1,---). 
Our goal is to examine the evolution of the system towards the steady-state by stopping 
calculations at the time derivative ed/dgi. Note that this mathematical ruse is widely used 
in literature, especially for solving the kinetic equations, expressed in natural guiding centre 
variables, valid within the drift approximation (see , for example t 



A quite similar procedure can be found in Ref. [15| for solving the drift-kinetic equation, /, 
for the toroidal rotation induced by neutral beam injection. In this case, the distribution 
function for fast ions is expanded in powers of the small parameter entering in the problem 
(the ratio of the slowing down rate to the bounce frequency) and the "singularity" 1/e is 
contained in the term at the lowest order, /_i, and not in the time derivative. Also in this 

12 



approach, the differential equation for the "reference state" (i.e. for /_i) is determined to 
the dominant order. The evolution on the various time-scales is obtained by applying the 
bounce averaging to the dynamical equations for each order in the perturbation expansion. 
With this method, the authors have been able "to capture" the main features of the slowly 
varying functions. The macroscopic observables, like the torque on the plasma, have then 
been obtained by applying the magnetic surface average. 

Coming back to our problem, condition i) is fulfilled if, during the relaxation, Eq. flT9|) is 
expanded starting from the term 

r« ^-1 = - (29) 



Taking into account Eq. (l25f) and Eq. ( 1271) . at the lowest orders we find 



dx;^" 



dg_ 



-6^0^'iVoA. [gS'J^. - x-^ 



(30) 



dXf^ dXp"^ 



i^A; 



, dg 



dgo 



nt^ AT I nl^^ T<^ V<^'^ 



with N^xn = M;i^6r,^ and Ml is the inverse matrix of M^^^ = 6i;+g^^X^-gor,.,x i.e., MlM^^ = 
6^. We note that Eq. (!30|) transforms in a covariant way under the following thermodynamic 
coordinate transformations (TCT) [l| 



X 



j: 



a/j. 



-j: 



X' 



Indeed, by inserting transformations ( 1311) into the Eqs (1301) we find 



(31) 



dQ-\ 



-eg't{9'o^''j'oZ-X'r)No,p 



dX'°'^ dX'°"^ 



(32) 



dX^ I dXf^ Jt^^l T.lni' fa v'aKX AT dX' 



ax"'' dx"!^ 



aA PiV' CiK 



dX'' 



or 



dX, 



l/iX AT/ / ~/Ki/ j'a 



l^"*^ dXf"'^ _ /^Atxt/ I ;^/ku j'a v'aK] 

\^ ds-i + deo ~ 9o ^^oxk\9o -^av ^o J 



(33) 

(34) 
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where 

Kxk = ^^"H^dx^dX^ ^^^^ 

Hence, matrix A^oak transforms like a thermodynamic tensor of second rank. The TCT 
are the most general forces-transformations leaving invariant the expression of the entropy 
production and the Glansdorff-Prigogine dissipative quantity V defined as 

dX"/^ 
V.J',— (36) 



Notice that from Eq.( l3T|) . the following important identities are easily derived 



Q 



yaL' "^ _ g . X''^'^ ——- = ("371 

aX'^'^dX'''' ' dX"'^dX"'^ ^ ' 

So the linear transformations 

^,a^, ^ cZ^x""" (38) 

where d^^ are constant coefficients (i.e., independent of the thermodynamic forces) are a 
special, but very important, class of the TCT. The covariance under translation of the 
variables X°"^ has been widely used in literature. We cite, as an example, the work of 



Hinton and Hazeltine IGj on the interpolation-expression (i.e., the expression valid for the 
fully collisional as well as the low-coUisional transport regime) for the diffusion coefficient 
of magnetically confined collisional-plasmas. The original interpolation formula, found by 
the authors by solving the gyrokinetic equations, showed a wrong sign as one goes from the 
low-collisional regime to the fully collisional regime. Hinton and Hazeltine have been able to 
restore the correct sign by performing a linear transformation of the thermodynamic forces. 
The linear transformations constitute a very important class of the TCT and it is always 
possible to associate to them a clear physical interpretation. For example, in the case of 
the Hinton-Hazeltine interpolation formula, the transformed flux corresponds to the total 
averaged radial energy flux obtained by the linear combination of the two original fluxes, the 
heat flux and the enthalpy carried out by the particle flux. However as already mentioned, 
the balance equations, as well as the closure equations, are covariant under the entire class 
of the TCT and, a priori, there are no reasons for retaining only the linear transformations. 
V can be obtained by deriving both sides of the flrst equation of Eqs ( l30l) with respect to 
parameter g^i. However, in accordance with the request ii) mentioned in section [Til this 
operation should be performed in such a way to preserve the covariance under TCT. For this, 
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we undertake the covariant differentiation along a curve of botli sides of the first equation 
of Eqs ([30]) 111 



2 "V^a/i 



d'X, 



dX^^ dX^'^ 



dgti 



dg^i dg-i 



e' U'^Jou - xn (go''' Jo.' - xtjgo'No,.' 



(39) 



X 



9'o"Nor,u]^ + nx.9o"No,' 



where comma (, ) stands for partial differentiation with respect to the thermodynamic forces 
and Tq^^ is the afiine connection evaluated at the lowest order. The elements F^^ transform 
like an afiine connection under the coordinate transformations Eqs fl3T|) . Its expression is 
reported in Ref. pj 



pM ^ }_„f^u f ^9\y dg^^ _ dgxK \ _}_^a^ij^au ^9>^ 



From Eq. f l25p . we obtain 



W-^X, ,_,) ^ ^ + F^,^^ ~ 0{e' 



dg_i 



dg_i dg^i 



This can be easily seen by noting that 



(40) 



(41) 



W"^(X, g_,) = W^{Xo, g-i) + e-H"^(Xo, X,, g_{) + 0{^ 



(42) 



where 



•H"'^(Xo,Xi,f?_i) 



rf^Xf^ , -p/. fdX^^dX^^ dX'^^dX;^^ 



pM 



+ 



(43) 



dg_i '^\dg^i dg_i dg^i dg_i 

From Eqs (l30|l we get 

?^"'^(Xo,Xi,f?_i)~0(e) (44) 

Thus, taking into account Eq. fl5U]) . close to the steady-state, the evolutionary-trajectory 
(ev.-tr.) of the system satisfies, at the dominant order, the equation 



W"^(Xe.._t..,f?_i) = + O(e' 



(45) 



We can check that Eq. (H5ll satisfies condition ii). Indeed, multiplying this equation with 
the thermodynamic forces X"'^ and contracting, we obtain 

£xy 



ja ^^ev.—tr. , Y°'l^ n 



jvau J van 

"■^ev.-tr. "'^ev.-tr. 

dg^-i dg^i 



+ O(e' 



(46) 
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Recalling that g^u — ^ 
following identity 



* 



n 

1, where ^ indicates the arc-parameter [1], we get the 



j: 



, " ^^ev.-tr. 



dP 



dq 



^^ev.—tr.y^iv,K, 



"'^^ev.-tr. ^^^ev.-tr. 



dxT 



(47) 



where P = J" j'"' *''• . Inserting Eq. P7|) into Eq. fj46l) . we obtain, at the leading order, the 
expression 



dP 



di, \2 



Oie' 



dg^i Kdg^i- 
Integrating from the initial condition to the steady-state, we get 

d^ ^2 



^\^ St. state) -* 



At the steady-state we have g_iP{X', 



a \ 

St. state/ 



V = g_^P = g_^J\ 



a ^ ev.—tr. 



, , , dg'_, > 

0, so we finally obtain 
d<; n2 



(48) 



(49) 



-Q-i 



dg'_i < 



(50) 



^ dg-i J \dg'_i- 

where the inequality is only saturated at the steady-state. Thus, after fast transient, up to 
e^, the correct (perturbed) evolution equations close to the steady-state are 






(fiX^^^ , in^ dX^^ dX^'' 



"dZ"i 



dXf^ 



Ak dQ—\ dg-i 



dXZ _^A«A; 



(51) 



nP-^ AT I T.K.i' ja vai 

9o ^^oak 9o -'ou - ^0 



^ dg^i dgo 

Eqs. 0511) are covariant under the TCT, satisfy the UCE and are derived from the balance 
equations Eqs (fT9|) . This set of equations are closedhj recalling the definition of the reference 
state. Hence, dX^^ /dg^ = 0. So, near the steady-states, the equations describing the 
relaxation of magnetically confined plasmas read (see also the footnote 23l |) 



( d'^X"" p^ dXg^ dXg" 
dg^ + -*- Ak dg dg 







(52) 



dX° 



/iX 



21 ^A*-^ AT I p.fii' 7Q VOiK 



where, for simplicity, we have suppressed in the parameter g_i the sub-index " — 1". Sys- 
tem ( l52l) admits a unique solution satisfying the initial conditions X^'^(O) = X"'^(0) [i.e., 
XnO) = 0] and XZ,state + eX^' 



st. state st. state' 
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It is known that, by a suitable choice of the parameter q, the differential equation for the 
shortest path in the thermodynamic space reads (see, for example, Ref. (l8| ) 

Tjaa(v ^^ — s.path -p/x s.path s.path ^ /rQ\ 

W iX,-patH,^) = ^^ + T,^^^ ^^-0 (53) 

Hence, from Eqs 0451) we straightforwardly obtain the following result: In the collisional 
transport regimes and under the validity of the drift parameter approximation, after fast 
transients, tokamak-plasmas relax to stable steady-states in the generalized frictions space so 
that 

J"W"'^(Xe.._i.., g) - J^U'^^iXs.path, q) = O(e') 

(54) 
(no summation over a) 

It is worthwhile mentioning that Eq. ( 154|) can be derived from Eq. (145|) . but the reverse is 
not true. Indeed, in our example we made use of the drift parameter approximation. In this 
particular case, parameter e is independent of the distance of the system to the steady-state 
and we obtained the result that, at the first order of approximation, collisional plasmas 
relax to the stationary states along the shortest path traced out in the space of the parallel 
generalized frictions. However, this result is not satisfied in general. The small parameter 
e, entering in a problem of relaxation depends , in general, on the distance of the system 



to the steady-state. It is possible to show that Eq. (I54p remains valid a 



so when e depends 



on the distance to the steady-state. In the paper reported in Ref. [llj, several examples 
of relaxation of chemical reactions to stable steady-states in the chemical affinities space 
are analyzed. Even in these cases, it is possible to show the validity of Eq. (^^ when the 
chemical systems relax near the stationary states. In these examples, X^ and J^ are the 
chemical affinities and the chemical reaction rates, respectively and e measures the distance 
of the system from the stationary state. 

Note that in Onsager's region, F^'^ -^ and (5'o'^A^oki/),a -^ 0. Moreover, ^ oc <j (see Ref. 
l|). Hence, the evolution equation for X^^ reduces to [see Eq. (139!) ] 

-^ = 0{e') or Xo"^(0=Xo"':,,,,,,+ (l-^)[x"'^(0)-Xo"':,,,„J+0(e2) (55) 

In other words, in the Onsager region, at the dominant order, the solution is a straight line 
in the space spanned by the parallel generalized frictions. 
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B. Spatially Dependent Magnetically Confined Plasmas 

In this case, firstly we have to develop the space-time dependent thermodynamic forces 
(rY"^), flows (JT"") and transport coefficients in (spatial) Fourier's series. Then, we perform 

ihe same calculations as in the homogeneous case by taking into account the slaving principle 

ITj l (see the footnote 2J|). Close to the steady-state, we have pj 






j;iv)it) = ^J^J^{r,t)exp{-zk.r)dv 

^JoW = l^fx"^ir,t)exp{-zk'-r)dv (56) 

dv denotes a (spatial) volume element of the system, and the integration is over the en- 
tire space V occupied by the system in question. The evolution equation for </(k)u ^^ ^^® 
thermodynamic mode (i.e., the mode with wave-number k = 0), reads 

rj^m, = -5,. {T^o)Jmu - X^o)) + O.T.(k, kOk^kVo + 0(6^) (57) 



■KU 
) 



where O.T.(k, k')k^k/^o stands for other contributions from different wave-numbers and 'g'^Q 
is a (time-dependent) positive definite matrix [1]. Due to the slaving principle, near the 
steady-state, the evolution equations for J/q> reduce to 

r J"(o)^ ^ -5m«(^(o)^°o). - X^o)) + O(e^) (58) 

By performing now the same multiple time-scale calculations as for the homogeneous case, 
we arrive to the following final result: the thermodynamic mode (k = 0) relaxes to the 
steady-state according to the law 

J("o),W"'^(X(o), g) - J^o)MZathiXio), q) = 0{e') (59) 



III. CONCLUSIONS 

In this concluding section, we shall use the results previously obtained for finally deriving 
the expression of the afiine connection. The Universal Criterion of Evolution theorem has 
been demonstrated by Glansdorff and Prigogine through the balance equations. From this 
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theorem and by analyzing examples of relaxation of chemical systems close to steady-states, 
we may obtain the expression for the afiine connection of the thermodynamical forces space. 
This construction may be made "step by step". We multiply the first equation of Eqs ( 15U]) 
by g^u and sum over the index fi. By taking the derivative, with respect to parameter g 
{g = Q-i), of both sides of the resulting equation, we find (for simplicity, we omit the index 
a) 

<PX>' 1 , . dX^ dX"" 

dg'^ 1 dg dg fQQ\ 

= e'{9oJl - X^') (go'^'Jow - Xf)g^mo,^,No.,x 
Thus, close to the steady-state, we have 

d^X" I, .dX^dX^ ^, 2N 

^-^^ + 2 (^^^- + ^-^)^^ = ^^' ^ ('') 

At the lowest order, previous equation may be rewritten as 

d^X" 1 „^, .dX^dX^ ^,2, . . 

^^+2*'"'(w+9..,.)^^ = 0(a (62) 

It is not difficult to show that Eq. (1621) satisfies the UCE but it is not covariant under the 
TCT. This because the term 1 /2g'^'^ (^g^Xf^ + g,yi^x) , does not transform as an affine connection 
under the TCT. This condition is satisfied by adding the Levi-Civita term —l/2g^'^gxK,u- 
Now, if one wants the Universal Criterion of Evolution satisfied also when the system relaxes 
along a shortest path, without imposing a priori any restrictions on transport coefficients, 
an extra term to the Levi-Civita affine connection should be added. It can be checked that 
the most general expression for this extra term is l/{2a)X^X'^gxK,iy [l]- Hence, the affine 
connection my be written as 

fM _ ^ u,y f 9guK dgxu dgxf, \ , 1 ^f^^^dgx^ 

^ ^'^ " 2^ Ux^ + dX- dX-J + 2a^ ^ dX- ^^^^ 

However, from the general theory on non-Riemannian geometry, we know that both equa- 
tions 

d'^X^' ~„ dX^dX" ^ .^, 

dg^ ^^ dg dg ^ ^ 

and 

d^Xt" „ dX^dX"" 

-nr + r^-^ — t- = o ^i^h 

dp^ dp dp 

TL = TL + S'x^^ + Ki^x (65) 
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are evolution equations for the same shortest path [iS!]. In Eq. ( 165|) . p is a parameter, related 
to g by 

— = cexpf — 2 / ipydX^] ; c = positive const. (66) 

and ipu is an arbitrary vector under TCT. In literature, modifications of the connection 
similar to Eqs (!65l) is referred to as projective transformations of the connection and ipK 
as the projective covariant vector. Therefore, the UCE is satisfied for every shortest path 
constructed with affine connections F^^, linked to F^^ by projective transformations Eqs (!65|) 
1181 1 ■ By standard methods in non-Riemannian geometry, it is not difficult to determine the 

n n 

final expression of the projective-invariant affine connection [1[, [18l |: 

^^ 2^ \dX- dX^ dX'^J 2a dX'^ 2{n + l)a[ ^ dX^ '^ dX'^ 



where n denotes the dimension of the generalized frictions space. We can check that the 
expression reported in Eq. fl67j) transforms under TCT (therefore, also under the translation 
of the variables X"^) as an affine connection. A non-Riemannian geometry has been succes- 
sively constructed out of the components of the affine connections [Ij. The main conclusion 
of the present analysis is thus: 

Close to steady-states, the geometry of the thermodynamic space is non-Riemannian with 
affine connection given by Eq. [6l\ ). The knowledge of the expression for the affine connec- 
tion is a fundamental prerequisite for the construction of the (nonlinear) closure theory on 



Q 



transport processes. In Ref . [l| , the curvature tensor and the nonlinear transport equations 
have successively been derived from Eq. fl671) and by introducing the following assumption: 
There exists a thermodynamic action, scalar under TCT, which is stationary with respect 
to arbitrary variations in the transport coefficients and the affine connection. 
From this principle, a set of closure equations, constraints, and boundary conditions have 
been derived. These equations determine the nonlinear corrections to the linear (" Onsager" ) 
transport coefficients. The nonlinear transport equations have successively been used for 
computing the particle and energy losses in magnetically confined plasmas and for deriving 
;he distribution density functions for the species a in several collisional transport regimes 
81] and Qj. We already mentioned in the introduction that, the thermodynamic forces are 
the variables, which characterize the departure of the system from the (thermodynamic) 
equilibrium. As a consequence, for plasmas not sufficiently close to the equilibrium, the 
distribution function for the fluctuations of the thermodynamic quantities deviates from a 
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Maxwellian and, in general, the closure equations may not be linearized with respect to 
the thermodynamic forces. From the microscopic point of view, this aspect can be easily 
understood in the following way. Let us indicate with {Y^} a set of n extensive variables 
describing the state of the plasma. We shall write the distribution function f{Y) for which 
{F^} lies between {Y^} and {y^} + {dY^} as (see also the footnote |25l |) 



f{Y)dY =f{Yi,Y2r--,Yn)dYidY2---dYn 

. r h^'^{Y)[Y,-Y,][Y,-Y,] ^^^^^^^ 
= cexp< — — — \dY1dY2 ■ ■ -dYn 



(68) 



where c is the normalization constant, Kb Boltzmann's constant and h^^^ the elements of a 
symmetric positive matrix, respectively. Y^ denote the mean values of Y^ 



% = / /(Fi, F2, ■ ■ ■ , Yn)Y^dY,dY2 ---dYn 



(69) 



Hence, the fluctuations ^j defined as 



^fi = Yf,- Y^ with (/i = 1, ■ ■ ■ , n) 



(70) 



obey the distribution function 

/(6,6,---^n) = cexp 



h^^'iOU, 



2K 



B 



(71) 



The following quantities will now be introduced 



X^'^K 



dlogf 



B- 



(72) 



From this definition we find 



2 <9^M 



(73) 



The variables X'^, referred to as the thermodynamic forces, clearly characterize the departure 
of the system from the equilibrium (notice that X^ — )■ as .^^ — ?■ 0, for yU = 1, ■ ■ ■ ra). Near 
equilibrium, the matrix h^^ is close to a constant (symmetric positive) matrix, say /ig'^, and 
the distribution function, Eq. f l68|) . takes the form of a Maxwellian. Thus, we can write 



X^ 



-K"^. 



/(O =icexp 



h^ f f 
2Kn 



(74) 



21 



for thermodynamic systems close to equilibrium. The conjugate thermodynamic flow are 
defined as 



■J.^-^ (75) 



dt 
We shall now assume that decay of a fluctuation follows the ordinary macroscopic law (see, 



for example, Ref. ll|) and write 

-Km. (76) 



dt '"^' 



a 



Near equilibrium, matrix fc^ is (almost) a constant matrix |^, say k^^, and Eq. fl75]) reduces 
to 

—jr = Jfj. — —l^Qpiiu (77) 

By combining now Eq. (1771) with the flrst equation of Eqs (17^ . we flnally obtain 

J^l - Lf.^X" with L^y = k^l^hour, (78) 

where /iqa*!/ denotes the inverse of the matrix /iq'^, i.e., hQ^ho^u = S^. In literature, matrix 
Lf^,y is called the Onsager matrix of the transport coefficients. The main conclusion of the 
present analysis is that: the transport relations can be linearized with respect to the ther- 
modynamic forces only when the system is sufficiently close to the (local) equilibrium (i.e., 
when the thermodynamic forces are "weak"). In this situation, the distribution function for 
fluctuations is a Maxwellian. Reversely, the transport relations are necessarily nonlinear in 
case of non- Maxwellian fluctuation- distribution function. Let us consider, for example, the 
case of magnetically confined tokamak-plasmas. After a short transition time, the plasma 
remains close to (but, it is not in) a state of local equilibrium. Indeed, starting from an 
arbitrary initial state, the collisions would tend, if they were alone, to bring the system 
very quickly to a local equilibrium state, described by a Maxwellian distribution function. 
But slow processes, i.e. free-fiow and electromagnetic processes, prevent the plasma from 
reaching this state. The distribution function for the fiuctuations of the thermodynamic 
quantities also deviates from a Maxwellian preventing the thermodynamic fiuxes from being 
linearly connected with the conjugate forces. 

We would like to end this section by adding some physical remarks. In this paper, we have 
shown that magnetically confined plasmas tend to relax towards the mechanical equilibrium 
following the shortest path, traced out in the space of the parallel generalized frictions. We 
may ask the following question: " why should the thermodynamic forces, and not other quan- 
tities, obey this law ? ". The answer is, because these quantities are purely non-conservative. 
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As known, the thermodynamic variables can be classified as conservative oi n on- conservative. 
A fiuctuation of a conservative variable can be "dissipated" only through the boundaries 
and, due to his severe constraint, its evolutionary-trajectory towards the steady-state may 
be very complex in the phase-space. On the contrary, a fiuctuation of a non- conservative 
variable can be dissipated freely into the surrounding and its evolutionary-trajectory tends 
to approximate that of the shortest path. We know that, in thermodynamics, one of the 
most dissipative quantity is V, the Glansdorff-Prigigine quantity, given by Eq. fj50|) 



where, we recall, q is the arc-parameter in the space of the parallel generalized frictions 
and g parameterizes the equation for the shortest path in this space. By noticing that 
the transport coefficients depend only on the thermodynamic forces X^, it turns out that 
also quantity V depend only on the parallel generalized frictions. Since the thermodynamic 
forces are non- conservative variables and quantity V is defined in the space of the parallel 
generalized frictions, it is not so surprising to have observed that, in the thermodynamical 
forces-space, the X'^ tend to follow the shortest path for reaching the steady-state. This law 
is so well satisfied experimentally that we can write 

j; = (7^.(X)X-'^ ^ ^;..(X)X"'^ [or, g^^iX) ^ ~g^,{X)] (80) 

during the gradual relaxation of the thermodynamic forces towards the (global) mechanical 
equilibrium, with 

■^^ St. state ~ 9tJ.u\Xst.state)Xg^^^^-f.^ = g ^^ [Xst. state jX^^^^^^^ [OT, g^u[Xst. state) = gfiu[Xst.state)\ 

(81) 
rigorously valid at the steady-state. We draw the attention to the reader that we prefer to use 
the expression "the system tends to follow the shortest path" instead of "the system follows 
the shortest path" because, as known, the Universal Criterion of Evolution, expressed by 
inequality fl79|) . cannot be derived from a variational principle [Ref. |l|, |2| and, our results 
Eqs ( l52l) ]. For this reason, in general, the shortest path is not the trajectory that minimizes 
the Glansdorff-Prigogine quantity. 
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